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ABSTRACT 



Context. Ground-based simultaneous multiband transit observations allow an accurate system parameters determination and may lead 
to the detection and characterization of additional bodies via the transit timing variations (TTVs) method. 

Aims. We aimed to (i) characterize the heavily bloated WASP-4b hot Jupiter and its star by measuring system parameters and the 
dependence of the planetary radius as a function of four (Sloan g', r', i', z') wavelengths and (ii) search for TTVs. 
Methods. We recorded 987 images during three complete transits with the GROND instrument, mounted on the MPG/ESO-2.2m 
telescope at La Silla Observatory. Assuming a quadratic law for the stellar limb darkening we derive system parameters by fitting a 
composite transit light curve over all bandpasses simultaneously. To compute uncertainties of the fitted parameters, we employ the 
Bootstrap Monte Carlo Method. 

Results. The three central transit times are measured with precision down to 6 s. We find a planetary radius Rp = 1.413 ± 0.020i?jup, an 
orbital inclination / = 88.°57±0.45° and calculate a new ephemeris, a period P = 1 .33823 144+0.00000032 days and a reference transit 
epoch Tq = 2454697.798311 ± 0.000046 (BJD). Analysis of the new transit mid-times in combination with previous measurements 
shows no sign of a TTV signal greater than 20 s. We perform simplified numerical simulations to place upper-mass limits of a 
hypothetical perturber in the WASP-4b system. 

Key words, transiting planets 



1. Introduction 

More than 100 extrasolar planets have been detected to pass in 
front of the disk of its parent stars since the first transit observa- 
tions reported nearly twelve years ago (Charbonneau et al. 2000; 
Henry et al. 2000; Mazeh et al. 2000). These close-in planets or- 
bit their host stars with periods typically smaller than ~ 10 days, 
probably formed at greater orbital distances and later migrated 
inward governed by processes which are still under debate. The 
discovery of each transiting extrasolar planet is of great interest 
for planetary science as these objects provide a unique access to 
an accurate determination of radii (down to a few percent) and 
masses via transit photometry and radial velocity measurements 
of the host star (Henry et al. 2000; Charbonneau et al. 2000). 
Moreover they allow one to plot thermal maps of the planetary 
surface via infrared spectra (Richardson et al. 2007; Grillmair et 
al. 2007, Knutson et al 2007), determine the planetary temper- 
ature profiles and permit studies of the stellar spin-orbit align- 
ment (Queloz et al. 2000, Triaud et al. 2010, Winn et al. 2011, 
Johnson et a. 2011). All of these form a set of astrophysically 
precious parameters, which are critical for constraining the for- 
mation and evolution of these interesting objects. Furthermore, 
the time intervals between successive transits are strictly con- 
stant if the system consists of a planet moving on a circular 
orbit around the parent star. If a third body is present in the 
system, it would perturb the transiting planet causing the time 



* Based on observations collected with the Gamma Ray Burst 
Optical and Near-Infrared Detector (GROND) at the MPG/ESO-2.2m 
telescope at La Silla Observatory, Chile. Programme 083.A-9010. 



interval between successive transits to vary. The resulting tran- 
sit timing variations (TTVs) allow the determination of the or- 
bital period and mass of the perturber down to sub-Earth masses 
(Miralda-Escude 2002, Holman & Murray 2005, Holman et al. 
2010, Lissaueret al. 2011). 

The transiting planet WASP-4b was discovered by Wilson 
et al. (2008) within the Wide Angle Search for Planets in the 
southern hemisphere (WASP-S, Pollacco et al. 2006). The planet 
is a 1.12 Mjup hot Jupiter orbiting a G7V star with a period 
of 1.34 day. WASP-4b was found to have a heavily irradiated 
atmosphere and inflated radius. Refined planetary orbital and 
physical parameters, based on transit photometry were presented 
by Gillon et al. (2009), who used the Very Large Telescope 
(VLT) to observe one transit, Winn et al. (2009), who mea- 
sured two transits with the Magellan (Baade) 6.5m telescope at 
Las Campanas Observatory, Southworth et al. (2009), who ob- 
served four transits using the 1.54m Danish Telescope at ESO 
La Silla Observatory and Sanchis-Ojeda et al. (2011), who ob- 
served four transits using the Magellan (Baade) 6.5m telescope. 
During two of the transits, Sanchis-Ojeda et al. (201 1) observed 
a short-lived, low-amplitude anomaly that the authors interpreted 
as the occultation of a starspot by the planet. Southworth et al. 
(2009) also noted similar anomalies in their light curves and 
the possibility that they were caused by starspot occultations. 
Sanchis-Ojeda et al. (201 1), combined their data set with that of 
Southworth et al. (2009) and found that each of them is consis- 
tent with a single spot and a star that is well-aligned with the 
orbit. Tracking this starspot, it was possible to measure the rota- 
tion period of the host star. A new possible starspot has been re- 
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cently detected in two closely spaced transits of WASP-4 by the 
MiNDSTEp collaboration with the Danish Telescope (Mancini 
2011, private communication). Furthermore, Beerer et al. (201 1) 
performed space-based secondary eclipse photometry, using the 
IRAC instrument on the Warm Spitzer Space Telescope, of the 
planet WASP-4b in the 3.6 and 4.5 yum bands. Their data sug- 
gest that WASP-4b's atmosphere lacks a strong thermal inver- 
sion on the day-side of the planet, which is an unexpected re- 
sult for an highly irradiated atmosphere. Caceres et al. (2011) 
analyzed high-cadence near-infrared ground-based photometry, 
detected the planet's thermal emission at 2.2 fim and concluded 
that WASP-4b shows inhomogeneous redistribution of heat from 
its day- to night-side. Finally, Triaud et al. (2010) investigated 
the spin-orbit alignment (J3) of the WASP-4 system by measur- 
ing the parent star radial velocity during a transit of its planet 
(Rossiter-McLaughlin effect) and found p - 4°^34„, excluding a 
projected retrograde orbit. 

In this paper we present three new transits of WASP-4b ob- 
served in August and October 2009 with the GROND instrument 
(Greiner et al. 2008), attached to the MPG/ESO 2.2m telescope 
at La Silla Observatory. We recorded each transit in four opti- 
cal g', r', i' and z' (Sloan) channels simultaneously. Using these 
new data sets we measure the planet orbital and physical param- 
eters and constrain its ephemeris by fitting a new photometric 
data set over four pass-bands simultaneously. Unlike all previ- 
ous observational studies related to the WASP-4b exoplanet, that 
fit data in a single bandpass, we model a transit light curve us- 
ing the four bandpasses simultaneously. As pointed out by Jha 
et al. (2000), this approach permits one to break a fundamental 
degeneracy in the shape of the transit light curve. Each transit 
light curve can be described primarily by its depth and duration. 
For a single band observations it is always possible to fit these 
with a larger planet if the stellar radius is also increased and if 
the orbital inclination is decreased. The main advantage of the 
multi-band transit photometry is that it allows the determination 
of the planet orbital inclination unique, independent of any as- 
sumptions about the stellar and planetary radii. Finally, we add 
the three new mid-transit times and search for TTVs. This pa- 
per is organized as follows. In section 2 we present an overview 
of the three transit observations of WASP-4b and the data reduc- 
tion. In section 3 we discuss the light curve analysis and the error 
estimation. Section 4 summarizes the main results. 



2. Observations and data reduction 

Three transits of WASP-4b were observed with the MPG/ESO 
2.2-m telescope at La Silla Observatory (Chile) during three runs 
on UT August 26 and 30 and October 8 2009 (see Fig. 1). Bad 
weather conditions prevented data collection during the remain- 
ing two nights allocated for the project. To monitor the flux of 
WASP-4, we used the Gamma Ray Burst Optical and Near- 
Infrared Detector (GROND). It is a gamma-ray burst follow- 
up instrument, which allows simultaneous photometric observa- 
tions in four optical (Sloan g', r', i', z') and three near-infrared 
i JHK) bandpasses (Greiner et al. 2008). In the optical channels 
the instrument is equipped with backside illuminated E2V CCDs 
(2048 X 2048, 13.5 jum). The field of view for each of the four 
optical channels is 5.4' x5. 4' with a pixel scale of 0".158 pixeL'. 

During each observing run we verified that WASP-4 and a 
nearby comparison star were located within GROND's field of 
view. The comparison star was ~ 1 .3' south-east from our target. 
For each run we obtained repeated integrations of WASP-4 and 
the comparison star for ~ 3.5 hr, bracketing the predicted central 



transit times. The exposure times were in the range 8-14 s, de- 
pending on the weather conditions (i.e. seeing and transparency). 
During each run we used the fast readout mode (~ 10 s), achiev- 
ing a cadence ranging from 18 to 24 s. To minimize system- 
atic effects in the photometry associated with pixel-to-pixel gain 
variations we kept the telescope pointing stable. However, due to 
technical issues caused by the guiding camera ~ 35% of the total 
number of images were acquired without guiding, that resulted 
in ~ 10-12 pixel drifts per observing run (~ 3.5 hr). 

In the Near-Infrared (NIR) channels we obtained images 
with integration time of 10 s. Longer exposures are not allowed 
due to technical constrains in the GROND instrument. It was 
found during the analysis that the JHK-images were of insuffi- 
cient signal-to-noise ratio (SNR) to be able to detect the transit 
and to allow accurate system parameter derivation from the re- 
sulting light curves. In addition, the pixel scale in the NIR de- 
tectors is 0.6 px"', which resulted in a poorly sampled star psfs 
(typically 4-5 pixels). As we carried out all observations in the 
NIR without dithering to increase the cadence and to improve 
the time sampling and the photometry precision, it is hardly pos- 
sible to stack groups of NIR images with prior background sub- 
traction. Therefore, we excluded the NIR data in the further anal- 
ysis. 

At the beginning of the observations on UT 2009 August 26 
and 30, WASP-4 was setting from an air-mass of 1.02 and 1.05, 
respectively which increased monotonically to 1.5 and 1.6 at the 
end of the runs. During the October 8 observations, the air-mass 
decreased from 1 .06 to 1 .02 and then increased to 1 . 15 at the end 
of the run. 

We employed standard IRAlQ procedures to perform bias 
and dark current subtraction as well as flat fielding. A median 
combined bias was computed using 22 zero-second exposure 
frames and 4 dark frames were used to compute the master dark. 
As non negligible fraction of the data was obtained without guid- 
ing, it was critical to perform the data reduction with flat fields 
of high quality. A master flatfield was calculated using the fol- 
lowing methodology. From a set of 12 dithered twilight flats we 
selected frames with median pixel counts in the range lOK to 
35K which is well-inside the linear regime of the GROND opti- 
cal detectors. Finally, the reduced (bias- and dark-corrected) flat 
frames were median combined to produce master flats for each 
of the g', r', i' and z' bands. 

Aperture photometry of WASP-4 and the reference star was 
performed on each calibrated image. To produce the differential 
light curve of WASP-4, its magnitude was subtracted from that 
of the comparison star. Fortunately, the later was of similar color 
and brightness, which reduced the effect of differential color ex- 
tinction. For example, the instrumental magnitude differences 
(comparison star minus WASP-4) measured on UT October 8 
2009 were ^g' = 0.587mag, ^r' = 0.669mag, A/' = 0.681mag, 
Az' = 0.687mag. As pointed out by Winn et al. (2009) the com- 
parison star is slightly bluer than WASP-4 (A(g' - /') - -0.094 
mag). 

To improve the light curves quality we experimented with 
various aperture sizes and sky areas, aiming to minimize the 
scatter of the out-of-transit (OOT) portions, measured by the 
magnitude root-mean-square (r.m.s.). Best results were obtained 
with aperture radii of 16, 20.5 and 17.5 pixels for UT 2009 
August 26 and 30 and October 8, respectively. 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Relative g', r' , i' and z' -band photometry of WASP-4 obtained with GROND. The transits occured (for each panel from top 
to bottom) on UT 2009 August 26, August 30 and October 8. The transit light curves obtained on each successive run are displayed 
with an offset of 0.025 in relative flux for a better illustration. 



The light curves contain smooth trends, most likely due to 
differential extinction. To improve its quality we decrease this 
systematic effect using the OOT portions of our light curves. We 
plot the magnitude vs. air-mass data for each run and chaimel 
and fit a linear model of the form: 

f{z)=a + bz, (1) 

where (a) and {b) are constants and (z) is the air-mass. Once 
the coefficients were derived, we applied the correction to both 
the transit and the OOT data. Performing this correction we no- 
tice that the baseline magnitudes of our target and the compar- 
ison star show sUght, nearly linear correlations with the {x,y) 
center positions of the PSF centroids on each CCD chip. We re- 
move these near-linear trends from WASP-4 and the reference 
star by modelling the base-line of the light curve with a poly- 
nomial that includes a linear dependency on the centroid center 
positions (x and y) of the stars on the detectors using the function 
of the form: 

f{x,y) = 1 + k_,x + kyy + k„xy, (2) 

where {kx), {ky) and {kxy) are constants. We then subtracted 
the fit from the total transit photometry obtained for each chan- 
nel during the three runs. No other significant trends that were 
correlated with instrumental parameters were found. 



Table 1. Light curve scatter difference measured with the stan- 
dard deviation of the time series prior and after detrending. All 
quantities in mmag. 



run 


g' 


r' 


i' 


z' 


Aug 26 


0.503 


0.358 


0.482 


0.339 


Aug 30 


0.818 


0.973 


0.633 


0.725 


Oct 10 


0.733 


0.614 


0.430 


0.495 



To illustrate the light curve quaUty improvement we present 
the scatter decrease after the described detrending procedure 
(Table 1). 

3. Light curve analysis 

To compute the relative flux of WASP-4 during transit as a func- 
tion of the projected separation of the planet, we employed the 
models of Mandel & Agol (2002), which in addition to the or- 
bital period P and the transit central time Tc are a function of five 
parameters, including the planet to star size ratio (Rp/Rt, where 
Rp and /?* are the absolute values of the planet and star radii), 
the orbital inclination i, the normaUzed semimajor-axis {ajR,, 
where a is the absolute value of the planet semimajor axis) and 
two limb darkening coefficients (mi and M2). Because the param- 
eters a//?„ Rp/Rt and i are degenerate in the transit light curve. 
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Fig. 2. The composite light curves of WASP-4 obtained during three runs on UT 2009 August 26, August 30 and October 8. The 
best-fit transit models are superimposed with continuous lines and the observed minus modeled residuals are shown centered at flux 
level 0.95 on each panel along a constant line. 



we assume fixed values for = 0.92 + O.O6M0 from Winn 
et al. (2009) and = 0.907;;};5[^/;. from Sanchis-Ojeda et al. 
(201 1) to break this degeneracy. Ideally we would like to fit for 
all parameters however, only the first three determine the best-fit 
radius of the star, the radius of the planet, its semimajor axis and 
inclination relative to the observer. 

The values for the two limb darkening coefficients were in- 
cluded in our fitting algorithm to compute the theoretical transit 
models of Mandel & Agol (2002). We assume the stellar limb 
darkening law to be quadratic. 



Table 

Law). 



2. Theoretical Limb-Darkening Coeflicients (Quadratic 



LD coefficient 


s' 


r' 




z' 


wavelength (nm) 


455 


627 


763 


893 


Ml (hnear) 


0.623 


0.413 


0.314 


0.248 


U2 (quadratic) 


0.183 


0.290 


0.303 


0.308 



To derive the best fit parameters we constructed a fitting 
statistic of the form: 



^ = 1 - Ml(l - yU) - M2(l - fJ-f, 



(3) = Z 



/•(observed) -/](predicted) 



(4) 



where / is the intensity and yu is the cosine of the angle be- 
tween the line of sight and the normal to the stellar surface. 
We fixed (initially) the limb darkening coefficients (mi and M2) 
to their theoretical values (see Table 2), which we obtained 
from the calculated and tabulated ATLAS models (Claret 2004). 
Specifically, we performed a linear interpolation for WASP-4 
steflar parameters = 5500 + 100 K, log g = 4.4813 + 0.0080 
cgs, [Fe/H] = -0.03 + 0.09 and Vt = 2.0 + 1.0 km s"', which we 
adopted from Winn et al. (2009). 



where /^(observed) is the flux of the star observed at the i-th 
moment (with the median of the OOT point normalized to unity), 
cr, controls the weights of the data points and /(computed) is 
the predicted value for the flux from the theoretical transit light 
curve. We assume the orbital eccentricity to be zero and mini- 
mize the statistic using the downhill simplex routine, as im- 
plemented in the IDL AMOEBA function (Press et al. 1992). 
The minimization method evaluates iteratively the statistic, 
yet avoiding derivative estimation until the function converges 
to the global minimum. 
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In the entire analysis we derive uncertainties for the fitted 
parameters using the bootstrap Monte Carlo method (Press et al. 
1992). It uses the original data sets (from each run and pass- 
band), with their A' data points to generate synthetic data sets 
also with the same number of points with replacements. We then 
fit each data set to derive parameters. The process is repeated 
until we get an approximately Gaussian distribution for each pa- 
rameter and take the standard deviation of each distribution as 
the error of the corresponding fitted parameter. 

Ground-based time series data are often proned with time- 
correlated noise (e.g. "red noise"; see Pont et al. 2006) and there- 
fore the data weights cr, need a special treatment, i.e. the uncer- 
tainties must be calculated accurately in order to obtain reliable 
estimates of the fitted parameters. It is a common practice to use 
the calculated Poisson noise, or the observed standard deviation 
of the out-of-transit data for the weights. Our experience shows 
that these methods often result in underestimated uncertainties 
of the modeled parameters. To estimate realistic parameter un- 
certainties we employ two methods. First we rescaled the pho- 
tometric weights o-j so that the best-fitting model for each band 
and run results in a reduced = 1, which therefore requires the 
initial values for the photometric uncertainties to be multiplied 
by the factors = l-396,Xred = 1-058 and;^^^^, = 1.179, for 
the UT August 26 and 30 and October 8 2009 runs, respectively. 

Second, we take into account the "red noise" in our data by 
following the "time-averaging" approach that was proposed by 
Pont et al. (2006) and used in the transit data analysis of various 
authors including Gillon et al. (2006), Winn et al. (2007, 2008, 
2009) and Gibson et al. (2008). The main idea of the method is to 
compute the standard deviation (scatter) of the unbinned residu- 
als between the observed and calculated fluxes, ctj and also the 
standard deviation of the time-averaged residuals, ctn, where the 
flux of each A' data points where averaged creating M bins. In 
the absence of red noise one would expect 



Table 4. Summary of the residual rm.s. (observed minus calcu- 
lated flux) of the unbinned and binned data. 



o-i / M 



(5) 



However, in reality ctn is larger than cri by some factor jS, 
which, as pointed out by Winn et al. (2007), is specific to each 
parameter of the fitted model. For simplicity we assume that the 
values of /3 are the same for all parameters and find its value by 
averaging over a range of bins with timescales consistent with 
the duration of the transit ingress or egress .i e. 10 - 30 min. The 
resulting values for jS are then used to rescale ctj in the;^f^ by this 
value (see Table 4) 

To derive transit mid-times {Tq) we fixed the orbital period 
P to an initial value taken from Sanchis-Ojeda et al. (201 1) and 
fitted the light curves from each run for the best Tq, planet to 
star size ratio, inclination and normalized semimajor axis. We 
then find the best orbital period P and Tq using the method dis- 
cussed in 4.2. Then we take the new values for the Tq and P and 
repeat the fit for the best planet to star size ratio, inclination and 
normalized semimajor axis. This procedure is iterated until we 
obtain a consistent solution. 



4. Results 

4. 1 . System parameters 

We set the mass and radius of the star equal to 0.92 ± 0.06 Mq 
and 0.907 ^Q ^}^ R^, respectively, and determine the best-fit radius 
for WASP-4b, Rp by minimizing the function over the four 
band-passes and the three runs simultaneously. Furthermore, 



run 


band 


cr 


>e 








[mag] 




[mag] 


Aug 26 


g' 


0.0025 


1.41 


0.0035 


Aug 26 


r' 


0.0025 


1.38 


0.0035 


Aug 26 


i' 


0.0028 


1.45 


0.0041 


Aug 26 


z' 


0.0028 


1.49 


0.0042 


Aug 30 


g' 


0.0023 


1.19 


0.0027 


Aug 30 


r' 


0.0017 


1.10 


0.0019 


Aug 30 


r 


0.0025 


1.23 


0.0031 


Aug 30 


z' 


0.0022 


1.29 


0.0028 


Oct 08 


g' 


0.0018 


1.22 


0.0022 


Oct 08 


r' 


0.0018 


1.17 


0.0021 


Oct 08 


i' 


0.0019 


1.29 


0.0025 


Oct 08 


z' 


0.0020 


1.33 


0.0027 



Notes, cr represents the residual scatter over the entire observing time 
interval on each run; cTsys is the rescaled value for the photometric 
weights, reflecting the presence of red noise. 



we assume that the inclination and the semimajor axis of the 
planetary orbit should not depend on the observed pass-band. 
Therefore, we constrained these parameters to a single univer- 
sal value. Following this procedure for the radius of the planet 
we found 1.413 ± 0.020 Rsup and the best-fit value for the or- 
bital inclination and semimajor axis / - 88.57° + 0.45° and 
a = 0.02300 ± 0.00036 AU, respectively. For comparison of 
the key parameter, Rp, Wilson et al. (2008) derived 1.416!° °^^ 
/?jup, Gillon et al. (2008) derived 1.304+[;{;54 ^^^^^ -^j^jj 
(2009) derived 1.365 ± 0.021 Rj^p, Southworth et al. (2009) de- 
rived 1.371+[; |Jg i^j^^ Sanchis-Ojeda et al. (2011) derived 
1.363 ± 0.020 ^jup. Our value for Rp is slightly higher than the 
cited radius in the literature and is dominated by the uncertainty 
of the stellar radius R,. The value for the orbital inclination, 
from our analysis is in a good agreement with earlier results. For 
example, Winn et al. (2009) found 88.56°;[J;^^, Sanchis-Ojeda 
et al. (2011) found 88.80°;^;^^, while Gillon et al. (2009) and 
Soutworth et al. (2009) found 89.35°^l% and 88.80° - 90.00°, 
respectively. The results for the fitted parameters of the transit 
hght curve modelling allow us to derive a set of physical proper- 
ties for the WASP-4 system. Table 3 exhibits a complete Ust of 
these quantities. 

To complete the light curve analysis we further fit the light 
curves allowing the linear and the quadratic limb darkening coef- 
ficients {u\ and M2) in each pass-band to be treated as free param- 
eters. We aim to compare the shapes of the transit light curves 
and the fitted system parameters from this fit and the curves de- 
rived using theoretically predicted limb darkening coefficients. 
In order to minimize the^^-statistic using all of the 1 1 parame- 
ters, including 8 limb darkening coefficients we again employed 
the downhill simplex algorithm and the bootstrap method. The 
new parameter values we find Rp - 1.409 ± 0.020/?jup, ; = 
88.75° + 0.40°, a = 0.02298 ± 0.00033 AU resuh in a slightly 
smaller value of the^^d function (1.001), as the fit is able to re- 
move some trends in the data. However, we report the parameter 
values derived using the theoretically predicted limb darkening 
coefficient as final results because the new fit is also more sensi- 
tive to systematics in the hght curve photometry. 

A set of multi-band transit hght curves also aUows one to 
search for a dependence of the planetary radius, Rp as a func- 
tion of the wavelength. Variations of Rp in some of the pass- 
bands could be produced by absorption hues such as water va- 
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Table 3. System parameters of the WASP-4 system derived from the light curve analysis. 



Parameter 




Value 


68.3 % Conf. Limits 


Unit 


Comment 


- Light curve fit - 












Planet-to-star size ratio 


P = Rp/R, 


0.15655 


0.00028 


- 


1 


Orbital inclination 


i 


88.57 


0.45 


deg 


1 


Normalized semimajor axis 


a/R. 


5.455 


0.031 


- 


1 


Transit impact parameter 


h = (a/R,) cos i 


0.136 


0.043 


- 


1 


Reduced chi-square 


2 


1.005 


— 


- 


1 


- Planetary parameters — 












Radius 




1.413 


0.020 


Rjup 


2 


Semimajor axis 


a 


0.02300 


0.00036 


AU 


2 


Mean density 


Pp 


0.582 


0.036 


Pjiip 


2 


Surface gravity 


gp 


15.36 


0.91 


m s"^ 


2 


Transit duration 


to 


2.1696 


0.0047 


hour 


1 


Ingress/egress 


'iiig/egr 


0.3014 


0.0039 


hour 


1 


- Stellar parameters - 












Mean density 


p. 


1.715 


0.029 


g cm"' 


1 


g' Linear limb-darkening coefficient 


Ml 


0.616 


0.046 




1 


g' Quadratic limb-darkening coefficient 


U2 


0.211 


0.060 




1 


r' Linear limb-darkening coefficient 


Ui 


0.427 


0.044 




1 


r' Quadratic limb-darkening coefficient 


U2 


0.303 


0.064 




1 


i' Linear limb-darkening coefficient 


Ui 


0.292 


0.047 




1 


i' Quadratic limb-darkening coefficient 


U2 


0.304 


0.060 




1 


z' Linear limb-darkening coefficient 


Ml 


0.241 


0.045 




1 


z' Quadratic limb-darkening coefficient 


M2 


0.386 


0.063 




1 



Notes. (1) Based on the analysis of the multi-band light curves presented in this work; (2) Values obtained using (1) and results for the M,, R,, 
Teff and Mp from Winn et al. (2009) and Sanchis-Ojeda et al. (201 1). 



Table 5. Best-fit radius in each band. 



band 


A 


Radius 




(nm) 


Rjup 


g' 


455 


1.409 ± 0.021 


r' 


627 


1.415 ±0.020 


i' 


763 


1.407 ±0.021 


z' 


893 


1.420 ±0.021 



Notes. The best-fit radii were derived after fixing the orbital inclination 
and normalized semimajor axis to their best-fit values from Table 3. 

por, methane, etc. in the planetary atmosphere. As a check for 
this assumption we fitted the data originating from each pass- 
band individually. Instead of fixing Rp/Rt, i and a/Rt to a a sin- 
gle universal value we set the planet to star radius ratio as a free 
parameter and keep the remaining quantities as free parameters. 
Table 5 displays the best-fit radii as a function of the wavelength 
with no indications of variations within the measured errorbars. 

4.2. Deriving the transit ephemeris 

One of our primary goals is to measure an accurate value for 
the transit ephemeris (Tq and P). We include all available light 
curves from the three runs and fit for the locations of minimum 
light using the best-fit planetary radius, inclination and semima- 
jor axis from Sect. (4.1.). We minimize the;^'^ as defined in Sect. 
(4) over the four-band data of each run by fitting simultaneously 
for the Tq. After the best-fit transit mid-times are derived we 
add them to all reported transit mid-times available at the time 
of writing in the literatur^ (Table 5). We further fit for the or- 

^ Southwort et al. (2009) measured mid-times during four transits of 
WASP-4b. We exclude these times in our analysis, as they were reported 



bital period P and the reference transit epoch To by plotting the 
transit mid-times as a function of the observed epoch {E) 

TciE) ^To + ExP (6) 

In this linear fit the constant coefficient is the best-fit refer- 
ence transit epoch (Tq) and the slope of the line is the planetary 
period, P. We then use the new values for the period and the 
reference epoch to repeat our fitting procedure for the system 
parameters until we arrive at a point of convergence. We used 
the on-line converter developed by Eastman et al. (2010) and 
transformed the transit mid-times from JD based on UTC to BJD 
based on the Barycentric Dynamical Time (TDB). We derive a 
Period = 1.33823144 ± 0.00000032 day and a reference tran- 
sit time Tc = 2454697.798311 ± 0.000046 BJD. The result for 
the best-fit period is consistent with the one from Sanchis-Ojeda 
et al. (2011), who found P = 1.33823187 + 0.00000025 day, 
and reported a smaller error as they measured the period using 
all available transit times, including the four additional measure- 
ments reported by Southworth et al. (2009). 

We further investigated the observed minus calculated (O - 
C) residuals of our data and the reported transit mid-times at 
the time of writing for any departures from the predicted values 
estimated using our ephemeris. As it was shown by Holman & 
Murray (2005) and Agol et al. (2005) a deviation in the (9 - C 
values can potentially reveal the presence of moons or additional 
planets in the system. We list and plot the 0-C values from our 
analysis in Table 5 and Figure 3, respectively. An analysis of the 
0-C values shows no sign of a TTV signal greater than 20 s, 
except the two big outliers at epochs 300 and 305. We put up- 
per constraints on the mass of an aditional perturbing planet in 



as unreliable due to technical issues associated with the computer clock 
at the time of observations (Southworth 201 1, private communication) 
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Table 6. Literature transit mid-times of WASP-4 and their resid- 
uals in addition to the ephemeris derived in this work. 



Epoch 



Transit midtime 
(BJD) 





300 
303 
305 
324 
549 
587 
809 
812 
815 
827 
830 
850 
859 



2453963 
2454364, 
2454368 
2454371 
2454396 
2454697 
2454748 
2455045 
2455049 
2455053 
2455069 
2455073 
2455100, 
2455112, 



5-7^5+0.0021 

<Qa4 1+0.00025 
9/:q 1-5+0.00097 

69623+°™" 

•"^"^■'-0.00026 

798228- 



O-C 

(day) 



Reference 



651228_QQQQQ72 

738643+o-«»'54 



.753274!! 
767816 

.826637_„;„„„„„, 
841 im+oiMSolo 

605y28_„ „„„„,. 
650002+''-''™ 

"-^^""^^-0 000074 



0.00011 
-0.00223 
0.00003 
-0.00171 
-0.00001 
-0.00008 
0.00012 
0.00016 
0.00010 
-0.00006 
-0.00001 
-0.00024 
-0.00005 
-0.00005 



Notes. (1) Gillon et al. (2009); (2) Winn et al. (2009); (3) Sanchis-Ojeda 
et al. (2011); (4) This work; The epoch of the first observed transit of 
WASP-4b was taken as equal to zero. 



the system as a function of its orbital parameters. We simplified 
the three-body problem by assumimg that the system is coplanar 
and initial orbits of both planets are circular. The orbital period 
of the perturber was varied in a range between 0.1 and 10 or- 
bital periods of WASP-4b. We generated 1000 synthetic O-C 
diagrams based on calculations done with the Mercury pack- 
age (Chambers 1999). We applied the Bulirsch-Stoer algorithm 
to integrate the equations of motion. Calculations covered 1150 
days, i.e. 860 periods of the transiting planet, that are covered by 
observations. The results of simulations are presented in Fig. 4. 
Our analysis allows us to exclude additional Earth-mass planets 
close to low-order period commensurabilities with WASP-4b. 

For the case of a transiting planet with a semi-major axis a\ 
and period Pi and a perturbing planet with semi-major axis 02 
on an outer orbit (i.e. ai > a\), period P2 and mass M2 Holman 
& Murray 2005 derived the approximate formula 



CHe = 



02(1 - ei) 



(7) 



(8) 



for the magnitude of the variation (in seconds) of the time 
interval (Af) between successive transits. One could imagine an 
exterior perturbing planet on a circular coplanar orbit twice as 
far as WASP-4b (period P ^ 3.75 day not in mean motion res- 
onance with the transiting planet). If such an imaginary planet 
had a mass of O.lMjup, it would have induced ~ 5 sec variations 
in the predicted transit mid-times of WASP-4b. Furthermore, the 
perturber could cause radial velocity variations of the parent star 
~ 13.78 m/s, which is below the radial velocity r.m.s. of the 
WASP-4 residual equal to 15.16 m/s, presented in the analysis 
of Triaud et al. (2010). Although, a few outliers are visible in the 
O-C diagram we consider the prediction of such an imaginary 
planet in the WASP-4 system as immature, because the weights 
of the O-C points might still require an additional term to ac- 
count for systematics of unknown origin. 



Oh' 

-3 'r 
-4L 







200 



400 
Epoch 



600 



800 



20 h 

10 h 

oh 



o -10^1 

° -2X>1 

-30 r 

-40 L 



550 



600 



650 



700 
Epoch 



750 



800 



850 



Fig. 3. The transit timing residuals for WASP-4b along with the 
one sigma errorbars. Top panel: The observed transit mid-times 
based on this work and others in the hterature were subtracted 
from the calculated times produced by our ephemeris. Lower 
panel: A closer view of the available transit mid-times from 2009 
and 2010. 




■ 10" 



- 10 



-= 10"- 



2.0 



5.0 



Pp^b 



Fig. 4. The upper-mass limit of a hypothetical additional planet 
that could perturb the orbital motion of WASP-4b as a function 
of ratio of orbital periods of transiting planet, Pb, and the per- 
turber. Pp. Orbits located in a grey area were found to be unsta- 
ble. 



5. Concluding remarks 

We have used the GROND multi-channel instrument to obtain 
four-band simultaneous light curves of the WASP-4 system dur- 
ing three transits (a total of 12 light curves) with the aim to re- 
fine the planet and star parameters and to search for transit tim- 
ing variations. We derived the final values for the planetary ra- 
dius Rp and the orbital inclination i by fixing the stellar radius 
Rt and mass M« to the independently derived values of Winn 
et al. (2009) and Sanchis-Ojeda et al. (2011). We include the 
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time-correlated "red-noise" in the photometric uncertainties us- 
ing the "time-averaging" methodology and and by rescaling the 
weights to produce value for the x^^i^ equal to unity. We further 
perform the light curves analysis by minimizing the function 
over all pass bands and runs simultaneously via two approaches. 
First we modeled the data using theoretically predicted limb- 
darkening coefficients for the quadratic law. Second, we fit the 
light curves for the limb-darkening. Both methods result in con- 
sistent system parameters within < 1 %. The second method pro- 
duced limb darkening parameters compatible with the theoreti- 
cal predictions within the one-sigma errorbars of the fitted pa- 
rameters. 

We added three new transit mid-times for WASP-4b, derived 
a new ephemeris and investigated the O - C diagram for out- 
lier points. We have not found compelling evidences for outhers 
that could be produced by the presence of a second planet in the 
system. We did not detect any short-lived photometric anomalies 
such as occultations of starspots by the planet, which where de- 
tected by Sanchis-Ojeda et al. (2011). At the transit r.m.s. level 
of our light curves (~ 2.2 mmag) it would be challenging to de- 
tect similar anomalies. However we note that due to the bright- 
ness of the parent star, the short planetary orbital period and the 
significant transit depth, the WASP-4 system is well suited for 
follow-up observations. 
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